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;  Abstract  . 

Using  the  inverse  scattering  transform  solution  for  parametric  amplification,  we 
show  that  the  result  of  the  interaction  of  a  small  signal  pulse  with  an  intense  pump 
wave  is  esronti^y  totally  independent  of  the  shape  of  the  pump.  Instead,  the  generated 
output  is  only  determined  by  the  amplitudes  of  the  TWI-sohtons  contained  in  the 
pump. 
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1  Introduction 


A  new  type  of  soliton  generated  by  a  quadratic  nonlinearity  has  recently  been  observed  in 
an  experiment  with  an  s3nichronously  pumped  optical  parametric  oscillator  [1].  Applications 
of  these  soUtons  for  second  harmonic  generation,  sum  frequency  generation  and  for  all- 
optical  switching  have  been  considered  in  references  [2]  -  [4].  In  Ref.  [3],  these  solitons 
were  termed  TWI-solitons.  An  analytical  theory  of  TWI-solitons,  based  on  the  inverse 
scattering  transform  (IST),  was  developed  more  then  20  years  ago  [5],  [6],  but  because  of  the 
mathematical  complexity  of  the  required  IST,  these  solitons  are  not  familiar  to  the  optical 
community.  In  this  letter,  we  use  the  TWI-soMton  theory  to  analyze  the  process  of  parametric 
amplification  of  ultrashort  pulses,  and  to  give  a  simple  interpretation  of  the  process. 

The  equations  describing  a  quadratic  interaction  of  three  waves  in  a  nonlinear  medium 
can  be  written  as: 
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where  Aj  are  the  normalized  ampUtudes:  Aj  =  {Ej/ Eo)^njUJ3jn^^,  ujj  are  the  frequencies, 
Vj  are  the  group  velocities,  Afc  =  —  ki  —  k2  is  the  phase  velocity  mismatch,  and  rij  are  the 

refractive  indexes.  Eo  is  a  normalization  amplitude,  defined  by  Eq  =  V (Stt )^, 
where  xm  is  the  coefficient  of  the  nonlinear  dielectric  susceptibility.  We  also  consider  the 
pulses  to  be  sufficiently  long  that  all  second-order  dispersion  effects  are  negligible  and  we  also 
take  the  ampUtudes,  Aj,  to  be  real.  In  this  paper  we  treat  the  case  where  the  group  velocity 
of  the  wave  with  the  highest  frequency.  U3,  lies  between  the  group  velocities  of  the  daughter 
waves,  vi  and  V2:  Vi  >  V3  >  V2.  This  case,  we  call  the  FSF  interaction  (Fundamental  - 
Sum  firequency  -  FVindamental)  interaction.  (In  Ref.  [7],  this  interaction  was  termed  the 
”  soliton-decay  case” ) . 

A  consideration  of  the  Eqs.  (1)  for  the  FSF  interaction,  based  on  IST  [5)- [8],  shows  that  if 
all  three  pulses  are  initially  separated  from  each  other,  then  imder  rather  general  conditions 
(such  as  the  ampUtude  never  crossing  zero.  etc.  [7]),  whenever  the  area  under  the  amplitude 
curve  of  any  of  these  waves  is  greater  than  7r/2,  then  that  wave  will  contain  solitons.  In  fact, 
the  number,  n,  of  soUtons  contained  in  that  envelope  follows  directly  from  the  area  of  the 
envelope: 

%  =  f  Aifl{t)dt  =  7rn  +  ei  (2) 

J -00 

where  Icj]  <  7r/2,  Aifi  is  the  amplitude  of  the  i-th  wave  and  n  is  the  number  of  TWI-solitons 
contained  in  the  wave.  The  quantity  €i  represents  the  non-soUton  (radiation)  part  of  the 
area. 

We  will  consider  only  the  case  of  parametric  amplification,  where  in  the  beginning  of 
the  interaction,  a  small  trigger  pulse  (a  signal  wave)  at  the  frequency  ui  is  placed  behind 
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some  intense  pump  pulse,  which  is  at  the  frequency  cJa  —  +  uii.  After  initialization,  the 

signal  wave  will  propagate  toward  the  pump,  will  overtake  it  and  then  interact  with  it.  Any 
time  the  high-frequency  pump  contains  a  soliton,  it  is  unstable.  When  it  decays,  it  emits 
exactly  one  soliton  into  each  daughter  wave,  for  each  soliton  in  the  pump.  Thus  each  soliton 
in  the  pump  can  be  thought  of  as  a  boimd  state  of  zero  binding  energy,  and  consisting  of  a 
signal  soliton  and  an  idler  soliton,  wherein  the  energies  of  each  are  perfectly  matched.  This 
matching  is  exactly  the  relation  that  is  prescribed  by  the  Manley-Rowe  relations.  When  this 
exact  balance  is  distmbed,  such  as  a  small  trigger  pulse  overtaking  the  pump  and  interacting 
with  it,  a  decoupling  of  the  signal  and  idler  solitons  occurs.  As  a  result,  the  pump  decays 
by  emitting  all  of  its  solitons  into  the  daughter  waves  (idler  and  signal)  [5],  [6]. 

As  it  follows  from  (2)  the  area  of  an  intense  pump  wave  (S3  >>  7r/2)  will  consist  almost 
exclusively  of  the  TWI-solitons.  Thus  we  shall  neglect  the  radiation  (non-soliton)  part  of 
the  pump  (since  ea  <<  S3),  and  from  now  on,  shall  assume  that  the  pump  wave  consists 
only  of  TWI-solitons. 

In  many  ways,  a  solution  of  nonlinear  equations  using  the  1ST  is  similar  to  a  solution  of 
the  wave  equation  iising  the  Fourier  Transform  [7|-[10].  In  order  to  solve  a  linear  differential 
equation  using  the  Fourier  transform,  one  has  to  decompose  an  arbitrary  shaped  wave  into 
a  superposition  of  simple  plane  waves.  The  amplitudes  of  such  plane  waves  do  not  change 
dvuing  the  propagation.  Therefore  the  spectrum  of  the  wave  at  the  end  of  propagation  is 
exactly  the  same  as  it  was  in  the  beginning.  All  that  has  changed  is  the  phases  between 
the  spectral  components.  To  solve  a  propagation  problem,  it  is  only  necessary  to  take  a 
superposition  of  the  original  plane  waves,  but  with  new  phases  as  determined  by  the  linear 
propagation. 

Analogous  to  the  linear  Fourier  transform,  in  the  nonlinear  1ST,  the  wave  with  an  ar¬ 
bitrary  shape  can  be  represented  as  a  nonlinear  superposition  of  sohtons  with  different 
amplitudes  and  phases.  Similar  to  the  linear  case,  the  number  of  solitons  does  not  change 
during  the  interaction,  but  their  phases  do.  Therefore  to  find  the  shapes  of  the  waves  after 
an  nonlinear  interaction,  one  needs  to  construct  the  final  envelopes,  using  the  same  solitons, 
but  with  the  new  phases  which  they  have  accrued  during  the  nonlinear  interaction. 

Let  A3  ^  be  a  time  profile  of  the  pump  before  the  interaction  (index  ‘i’  here  stands  for 
‘initial’).  Let  us  also  suppose  that  the  pump  initially  consist  of  n  TWI-solitons  with  the 
amplitudes  2r]i  and  the  phases  Di-  Because  of  the  nonlinear  nature  of  the  three- wave  inter¬ 
action,  a  nonlinear  superposition  can  not  be  obtained  merely  by  adding  together  all  solitons 
contained  in  the  pump.  For  the  TWl  interaction,  the  law  of  superposition  of  n  TWI-solitons 
is  given  by  what  is  called  the  n-soliton  formula.  For  TWI-solitons,  this  is  [7]: 

j,k=i 


where  /  is  the  identity  matrix  and  the  matrix  N  is  given  by 


Djexpl-(rjk+Vj)t] 

Vk+rjj 


(4) 
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In  (3)  the  negative  power  denotes  the  inverse  matrix,  and  2r]i  are  the  amplitudes  of  the 
soli  tons  contained  in  the  pump.  As  can  be  seen  from  (3)  and  (4),  the  residt  of  a  nonlinear 
superposition  of  n  TWI-solitons  depends  not  only  on  the  amplitudes  of  solitons  2T]i,  but  also 
on  their  phases,  Di.  (Strictly  speaking,  the  coefficients  Dj  are  not  the  phases,  but  they  are 
functions  of  the  phases.  Still,  we  shall  refer  to  them  loosely  as  the  ’’phases”.) 

When  the  TWI-solitons  are  far  apart  from  each  other,  the  coefficients  A  have  a  simple 
physical  meaning:  they  determine  the  positions  and  phases  of  the  TWI-solitons.  To  xmder- 
stand  this,  let  us  consider  a  simple  case  when  the  pump  contains  only  one  soliton.  For  ti=1, 
(3)  and  (4)  give: 

Ag  ^(t)  =  2T]isgn{Di)sech[2r)i{t  -  to)]  (5) 


where  to  =  ln(|I>i|/2j7i)/2r;i.  As  can  be  seen  from  (5),  to  gives  the  position  of  the  center 
of  the  soliton  with  respect  to  the  origin.  However  this  simple  physical  interpretation  of  the 
phases  A  is  legitimate  only  when  the  solitons  are  well-separated  from  each  other.  When 
the  soUtons  are  too  close  to  each  other  and  overlap,  then  it  does  not  make  sense  to  assign 
a  position  to  each  individual  soliton  within  the  group,  since  the  individual  solitons’  shapes 
merge  and  ”  mix” . 

Let  us  suppose  that  initially  there  were  only  the  pump  and  a  small  signal  wave,  and  that 
the  idler  wave  is  absent  before  the  interaction.  During  the  interaction,  all  solitons  in  the 
pump,  will  transform  into  idler  and  signal  wave  solitons  [5],  [6].  The  number  of  sohtons  in 
each  of  these  waves  will  therefore  be  equal  to  the  number  of  solitons  initially  contained  in 
the  pump.  However  the  phases  of  each  soliton  will  be  changed  by  the  interaction.  (If  the 
phases  did  not  change,  the  shapes  of  the  signal  and  idler  waves  on  the  output  would  exactly 
repeat  the  shape  of  the  pump  before  the  interaction).  The  final  phases  in  these  envelopes 
can  be  easily  obtained  using  the  analytical  theory  in  Ref.  [7].  For  the  initial  conditions 
considered  above,  after  some  algebra,  one  obtains  the  following  expression  for  the  phases 
D[^1  (the  index  ‘/’  stands  here  for  ‘final’)  of  the  k  -  th  soliton  in  the  signal  wave  after  the 


interaction: 


^(/)  ^  2pi{r]k)(nik 
iqn  Vk-om, 


exp 


(6) 


The  index  ‘1’  here  corresponds  to  the  signal  wave,  q  =  {vi—V3)/{vi  —  V2),  2r}k  is  the  amplitude 
of  the  k  -  th  soliton  in  the  signal  pulse  (which  is  a  times  the  corresponding  amplitude  in 
the  pump),  p{r)k)  is  the  reflection  coefficient  of  the  1ST  theory  for  the  initial  signal  wave, 
and  T]j  are  the  77’s  that  were  initially  in  the  pump  pulse.  When  the  intensity  of  the  signal 
pulse  is  small  initially,  then  we  have  p{7)k)  —  Aifi{t)exp{2T]kt)dt,  where  Ai,o  is  the  initial 
amplitude  of  the  small  signal  wave,  and  provided  the  integral  makes  sense. 

Formula  (6)  shows  a  remarkable  quality  of  the  soliton  interaction:  the  final  phases,  D\^k> 
depends  only  on  the  initial  soliton  content  of  the  pump  and  do  not  depend  on  the  initial 
phases  between  the  solitons  in  the  pump  ^3  *..  This  means  that  the  result  of  the  nonlinear 
interaction  does  not  depend  on  the  exact  pump  shape,  but  depends  only  on  the  soliton 
content  of  the  pump.  An  illustration  of  this  fact,  obtained  by  a  numerical  calculation  of  the 
system  (1),  is  given  in  Fig  1.  The  first  three  pictures  in  this  figure  show  the  shape  of  the  large 
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pump  and  the  small  signal  pulses  before  the  interaction.  Although  the  shapes  of  the  pump 
pulses  in  Fig  1  (a,b,  and  c)  are  very  different,  they  all  contain  exactly  two  TWI-solitons,  with 
r/’s,  7/1  =  1  and  r}2  =  2.7.  Because  of  the  identical  soHton  content,  all  three  waves  in  Fig  1 
(a,b,c)  have  the  same  energy  and  the  same  area  under  the  amphtude  curve.  The  expression 
(6)  predicts  that  the  output  for  the  signal  wave  will  be  exactly  the  same  for  all  three  cases. 
The  numerical  computations  are  in  complete  agreement  with  this  conclusion:  the  output  for 
all  three  processes  is  shown  in  the  Fig  1  (d),  which  consists  of  only  one  of  the  results,  with 
the  other  two  being  so  close  to  this  one,  that  no  differences  in  the  graphs  could  be  seen.  As 
can  be  seen  from  this  pictme,  the  two  solitons  which  initially  belonged  to  the  pump  finally 
moved  to  the  signal  and  idler  envelopes. 

Let  us  consider  the  following  physical  consequence  of  (6).  The  sign  of  this  expression 
is  determined  by  the  product  in  the  denominator.  Without  loss  of  generality  we  can  take 
'Hj  ^  Vj+i-  largest  soliton,  this  product  is  positive.  After  that,  will  change 

sign  every  time  j  is  increased  by  one.  Now,  according  to  (5),  this  means  that  the  largest 
TWI-soUton  in  the  signal  wave  will  have  a  positive  phase,  the  next  TWI-soliton  will  have  a 
negative  phase,  and  the  sign  will  change  from  one  soliton  to  another,  as  one  can  clearly  see 
in  the  example  of  the  numerical  calculations  for  two  TWI-solitons  in  Fig  1  (d). 

There  are  further  physical  consequences  of  (6) .  Consider  the  position  of  each  TWI-soliton 
in  the  output  pulse.  We  can  estimate  those  by  taking  the  logarithm  of  and  dividing  the 
result  by  2r]k'. 

(7) 

Vi  27?fc 

where  in  <5(j7fc),  we  have  collected  all  the  terms  which  depend  only  on  r]k  in  a  nonsingular 
fashion.  The  first  term  in  this  expression  shows  that  all  TWI-solitons  in  the  final  signal 
envelope  move  with  the  same  speed  v\.  The  second  term  in  (7)  shows  how  much  this  TWI- 
soliton  will  be  delayed  by  the  nonlinear  interaction,  this  term  will  be  negative  for  small  Ai,o- 
The  smaller  the  trigger  pulse,  the  longer  the  delay.  The  small  solitons  corresponding  to  the 
smaller  T)k  experience  the  larger  delays.  Therefore,  the  smaller  the  amplitude  of  the  signal 
wave,  the  farther  apart  from  each  other  will  be  TWI-sohtons  in  the  output  signal  wave. 

In  conclusion,  using  an  analogy  between  the  nonlinear  1ST  and  Hnear  Fourier  transform, 
we  represent  an  intense  pump  pulse  as  a  nonhnear  superposition  of  the  TWI-solitons.  We 
have  shown  that  the  result  of  the  process  of  parametric  amplification  of  the  small  signal 
wave  by  the  intense  pump  does  not  depend  on  the  shape  of  the  pump.  The  resulting  train  of 
TWl-soUtons  in  the  output  signal  wave  is  determined  only  by  the  number  and  the  intensity 
of  the  TWI-solitons  initially  contained  in  the  pump.  Simple  calculations  show  that  the 
separation  between  the  emerging  solitons  become  larger  when  the  intensity  of  the  trigger 
pulse  is  small.  The  numerical  calculations  are  in  a  remarkably  excellent  agreement  with  the 
analytical  theory. 
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3  Figure  captions 

1.  Figure  1.  Numerical  verification  of  analytical  calculations.  Pictures  (a,b,c)  represent  a 
temporal  distribution  of  the  pump  wave  calculated  using  (5).  Each  pulse  contains  exactly 
two  TWI-solitons  with  the  amplitudes  tji  =  2.7  and  772  =  1,  but  taken  with  different  phases: 
Di  =  2,  D2  =  1  (a),  =  5,  £>2  =  1  (b),  Di  =  -3,  £>2  =  1  (c).  The  output  of  the  nonhnear 

interaction  for  all  three  cases  is  given  in  the  Figl(d):  the  signal  wave  is  on  the  left  and  the 
idler  wave  is  on  the  right.  Both  TWI-solitons,  initially  contained  in  the  pump,  now  belong  to 
the  signal  and  idler  waves.  Since  the  pump  did  not  contain  any  radiation,  it  has  completely 
disappeared  by  the  end  of  the  interaction. 
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